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Abstract 

Implication of our recent proposal [J. Phys. Soc. Jpn. 65 (1996) 687] to 
treat large-amplitude gauge-field fluctuations around the slave-boson mean- 
field theory for the t-J model has been explored in detail. By attaching 
gauge flux to spinous and holons and then treating them as free gf-on's 
which respect the time-reversal symmetry, the optimum exclusion (g) and 
exchange {a) statistics have been determined in the plane of doping rate 
and temperature. Two different relations between a and g have been inves- 
tigated, namely g = \a\ (Casel) and g = \a\{2 — \a\) (Case2). The results 
indicate that slave fermion is favored at low doping while slave boson at 
high doping. For two dimension, in Casel intermediate statistics are found 
in between, while in Case2 no intermediate statistics are found. The conse- 
quences of varying the dimensionality and strength of J have been studied 
also. The latter has no qualitative effect for both cases, while the former 
has a profound effect in Casel. 

KEYWORDS: t-J model, gauge-field fluctuation, exclusion statistics, 

^-on 



1. Introduction 

Since the discovery of high-Tc cuprates-^^ there are many attempts to 
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understand this interesting class of materials. It has generally been ac- 
cepted that the essence of electronic properties is in the Cu02 layers, and 
that the low-lying excitations of the layers are described by a single-band t- 
J modeP'^^. Though the exact results are not known, the slave-boson mean 
field theory^"^) has played important roles in offering explicit predictions to 
be compared with experiment as regards the doping dependence and effect 
of fermi surfaces, and the agreement between the theoretical results and 
what have been seen in the experiments on the spin excitations probed by 
neutron scattering^^^ and nucler magnetic resonance (NMR)-'^-'^^ phonons^^^ 
and transport properties'"^-* are noteworthy. The scheme is based on spin- 
charge separation and the electron operator is decoupled into spinon and 
holon, which is assumed to be fermion and boson, respectively. The phase 
diagram derived by the mean field approximation in this scheme is given 
by Fig.l. The states realized include antiferromagnetic (AF) state, super- 
conducting (SC) state, "spin gap" (SG) phase, anomolous metallic (AM) 
state, and electron liquid (EL) state. The AM state, where both spinous 
and holons move coherently and essentially independent of each other, i.e. 
spin and charge are separated, is described as the uniform RVB (i/,-RVB) 
state. In the SG state, where the spin and charge are also separated, the 
spinous make short range singlet pairs, which is described as the singlet 
RVB (s-RVB) state. On the other hand in the EL state spinous and holons 
are no longer independent, because of the bose condensation of holons, and 
coupled together to form electrons. Based on this phase diagram various 
physical quantities especially, the spin gap phenomenon typically seen in the 
rate of NMR^'^ have been discussed^"^) in connection with the SG state in 
the phase diagram. 

It is generally the case, however, there exist fluctuations around the 
mean field solution. The most important fluctuations, in the present case, 
are expressed by the phases of the order parameters in the li-RVB state, 
which are viewed as U(l) gauge fields''*^. The charge transport properties 
governed by holons are essentially influenced by these gauge fields because 
of the small characteristic energy of holons. The linear T-dependence of 
the resistivity has been discussed by treating the gauge field in the Gaus- 
sian approximation'^^. However there are indications that the gauge-field 
fiuctuations can be large. For example, high temperature expansion on the 
lattice'^^ indicates that for a reasonable value oit/ J the root-mean-square 
of fiuctuations of gauge fiux is large (~ tt per plaquette) . This might lead 
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to quantitative^^) and even qualitative^^) changes in the mean field phase 
diagram. The existence of such large fluctuations may indicate the necessity 
to elaborate the description of the ground state, and there are some efforts 
towards this target. On this line, there are two approaches both of which 
try to incorporate the fact that in the limit of vanishing doping rate, i.e. 
half-fllling, the 7r-flux state is the most stable state. One idea exploits the 
equivalence of the 7r-flux state to d-wave due to SU(2) symmetry^^'^'^\ while 
another maps the 7r-flux state into slave semion scheme^ by binding fluxes 
of half-quanta to both spinous and holons. The latter type of approach 
to treat strongly correlated electron systems may date back to the work of 
Kalmeyer and Laughlin who mapped the two-dimensional spin system 
to the fractional quantum Hall system. However, since flux state is stable 
only at half-filling, it is an open question as regards what will happen in 
finite doping rate. 

Recently, we have proposed a different approach^"^) to treat these large 
gauge-field fiuctuations. The spinous and holons may avoid these large 
gauge-field fiuctuations by partly absorbing them, which results in statis- 
tical transmutation of these particles. The statistics of holons and spinous 
will be chosen such that the free energy is minimum, i.e. the statistics are 
to be optimized. Here the time-reversal symmetry is properly taken into 
account by considering the fiuctuating nature of the gauge field. A similar 
problem of determining variationally the optimum statistics of spinous and 
holons has been studied also by Mori^^) for the case of the Hubbard model 
with infinite Coulomb interaction. His formalism is, however, limited to 
continuum space, breaks time reversal symmetry and it turns out that only 
slave boson is realized at all doping rates. 

In our explicit calculations of the thermodynamic properties the ex- 
clusion properties of anyons (particles with flux) has been noted. This 
exclusion property is based on the idea of Haldane^^) which generalizes the 
concept of Pauli exclusion principle. A particle obeying this exclusion prin- 
ciple is called "^-on"^^) or "excluson" . Its distribution function was de- 
termined by Wu.^^) Recently, some models are known to realize ^-on's^^'^°). 

This paper has the following outline. Section 2 describes our model and 
formulation. The results of our numerical calculation are given in Sec. 3. 
In Sec. 4 conclusions are given together with discussions on the possible 
relevance in describing the anomalous metallic states of high-Tc cuprates. 
A part of the results for two dimensions in this study have been published 
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elsewhere^^^ . 

2. Model and Formulation 

2.1) Slave Boson Mean-Field Theory 

The t-J model is defined by the Hamiltonian, 

H = -Y^ [tij{l - ni_^)cl^Cj^{l - nj_^) + h.c] 

+J ^ SiSj - fie^rii^. (2.1) 

<i,j> i,a 

where the parameters and operators, especially, = \c\^Ta^Cip with r 
being Pauli's spin matrices, have the usual definitions, and <, > means the 
summation over pairs of nearest neighbor sites. 

In slave boson scheme the electron operator is decomposed in terms 
of holon {hi) and spinon (s^o-), c^fj — h^Sifj. The holons and spinons are 
assumed as bosons and fermions, respectively. The Hamiltonian is rewritten 
as 

^ = - E [tijhislsj^h] + h.c] +J ^i^j 

i,3,o- <i,3> 

-Me ^l^- - E (E ^l^- + - 1) (2.2) 

i,a i cr 

where the last term with Lagrange multiplier, A^, excludes double occupancy 
at any sites. In the mean field approximation the local constraint is being 
relaxed to global one, — > A. It has been shown^^ that the choice of values 
of the transfer integrals, {t/J = 4) t' /t = -|, t"/t = and t'/t = -|, 
t" /t = I simulate the fermi surfaces of LSCO and YBCO, respectively. 
Here t, t' , and t" are transfer integrals between the first nearest neighbors, 
the next nearest neighbors, and the 3rd nearest neighbors, respectively. 
For simplicity our work will consider only the standard t-J {t/J = 4 and 
t'/t = t" /t = 0) since the present conclusions are expected to be insensitive 
to these choices. Moreover, although the mean field theory predicts various 
states as shown in Fig.l we will limit ourselves to the w-RVB state, where 
Xs =< Xla Aa^i(^ ^ X/i =< hih^- > are the order parameters. The 
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effect of singlet spinon pairing will be studied in a separate publication. 
Then, the mean field hamiltonian, Hmf, is 

<i,j>,a i,a 

- E + h.c] -l^hJ2 ^l^i (2.3) 

<i,3> i 

where fih = \ ijls = l^e + IJ^h are the chemical potentials of holons and 
spinous, respectively, and Ujih) = txsi Uji^) = tXh + 3Jxs/8 are the effec- 
tive hopping integrals of holons and spinous. 

2.2) Fluctuations around the Mean-Field Solution 

The original slave boson scheme without mean field approximation the 
electron operator, Cia = /ijsio-, is invariant under local phase transformation 
which is lost totally in the mean field approximation. As Baskaran and 
Anderson-*^"^) noted, to recover this local symmetry, the fluctuations around 
the mean field solutions should be treated as U(l) gauge field and hence 
the order parameters are not constants but fluctuating in space and time; 
tij{q) = tqC^"'^^ (q = h, s) where tg are the mean field values. The fluctuation 
ttij is, indeed, a gauge field that affects the spinous and holons together. 
This cures the complete decoupling of the electron operator at the mean 
field level. The Hamiltonian will then be written as 

<i,j>,a i,a 

-^h ¥"''h\^3 + h-c] -^JihY, h\hi. (2.4) 

<i,j> i 

This is invariant under the following gauge transformation, 

aij — y ttij + (0i - 0j). (2.5) 

In the calculation of transport coefRcients^^^ the fluctuations have been 
treated in the Gaussian approximation and it has been indicated that gauge 
fluctuations are not small. 
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Here we propose a way to take account of these large gauge-field fluc- 
tuations by assuming that the fluctuating gauge fluxes are approximately 
of Ising type, i.e. the magnitude is nearly constant but its sign is changing 
spatially and temporally [Fig. 2(a)]. In this paper we will determine the 
magnitude of these Ising-type fluctuations that optimize the free energy by 
exploiting the ideas of flux binding and exclusion statistics. 



2.3) Flux Binding and Exclusion Statistics 



In a region where gauge-field fluctuations have some fixed value and 
then the uniform (internal) magnetic field is present, the fluxes due to 
gauge- field fluctuations will be attached to spinous and holons [Fig. 2(b)]. 
Now we have a system of anyons whose exchange statistics a (defined as 
the phase factor e'^^'^ accumulated by the total wave function when two 
particles are interchanged) has approximately a constant absolute value. 
The new particles are then defined as Sia = Siae'^^^'^ and hi = hie'^^^\ 
Here Qi = Y^j^iOi{i){h]hj + s]„Sja) and 9i{j) represents the angle of 
a vector from site j to i with respect to some flxed reference axis. It is 
straightforward to see 

hji] = e'^'^'hyhi, hihj = e'^'^hi, i ^ j 

and 

and then the exchange statistics of holons and spinous are a and ±1 + a, 
respectively. The onsite commutation relations are simply given by the 
commutation relations of the undressed particles. The Hamiltonian for hi 
and Sia- is the same as Eq.(2.4) but we expect that the magnitude of fluc- 
tuations dij = aij — a{Qi — Qj) are much reduced. Then we assume that 
the fluctuations dij felt by the new particles hi and Sia are negligibly small 
and the mean fleld approximation is meaningful. If we flx a, it will break 
the time-reversal symmetry. However, in the present case the time-reversal 
symmetry is maintained due to the fluctuating nature of a; the magnitude 
is flxed but the sign is fluctuating as shown in Fig.2(c), so that the (net) 
flux will vanish as an average and the thermodynamic properties of the 
system will be purely described by the exclusion properties by the dressed 
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holons and spinons. In other words, although the system is composed of 
domains with fluxes a and —a, it will be regarded as homogeneous in its 
exclusion properties if we neglect domain boundaries [Fig. 2(d)], and its ther- 
modynamics may be approximately treated by use of the ^r-on distribution 
function^®) 

ng{e) = {wg{e) + g)-'^. 
Here Wg{e) = w is determined by 

w^il^wy-^ = e^^'-''\ (2.7) 

with one-particle energy £, chemical potential fi and inverse temperature 
f3. For a special value of = (^f = 1) we have the familiar bose (fermi) 
distribution function. The distribution functions for several choices of g are 
shown in Fig. 3. Here the exclusion statistics is deflned as the decrease in 
the size of one-particle Hilbert space as the number of particles increases^^^ . 
For boson this size does not change, g = 0, and for fermion the decrease is 
one whenever the number of particles increases by one, g = 1- In between 
boson and fermion we have g-on, < g < 1. 

There are so far two proposed relations between the exclusion statistics 
g{>0) and exchange statistics a {—1 < a < 1); namely, 

Casel : = |ck| (2.8a) 

and 

Case2: g = \a\{2 - \a\) . (2.86) 

These relations, Eqs.(2.8a) and (2.8b), are the exclusion statistics of system 
of free-anyon with Dirac-type dispersion^^^ and with quadratic dispersion^^\| 
respectively, which are classifled as Casel and Case2 in the following. In 
the present study to see overall tendency as the doping rate is varied, the 
energy spectra of spinons and holons vary in a wide range and the correct 
assignment is not known and then we simply assume either (2.8a) or (2.8b) 
and investigate the consequences for both cases. If we assign the exclusion 
statistics gh = g{o!.) for holon, then gs = g{l — a) for spinon since they 
are composite particles of electron. Hence the free energy becomes a func- 
tion of 5 (doping rate), T (temperature) and a (statistics). By treating 
a(0<Q;<l)asa variational parameter we can minimize the free energy 
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and determine the phase diagram in the plane of 6 and T. In this study the 
dimensionahty, d, of the system can be arbitrary and we will examine the 
cases of = 1, 2, 3. 

The free energy is given by 

F = -2r^ln(l + - T^ln(l + w;^\u^)) 

k k 

3 J 

+iV [d{2tXsXh + yXs) + (1 - S)fis + Sfih] (2.9) 
and the self-consistent equations are 

and Xs = -nTY.,lk^gsi£k), (2-10) 



Xh 



where e-^ - -'^tgjj; - fJ^s, - -'^thlk ~ l^h with = Z]a;=i coskx, Ug^ 
{q — h, s) the "^r-on" distribution function and is the total number of 
lattice sites. 

3. Results 

In the following two subsections the dimensionality of the system is 
fixed to two, d = 2, whereas effects of the dimesionality will be discussed 
afterwards. 

3.1) T = 

First we will consider the T = where the physics and the overall 
features of the results are transparent and easy to comprehend. In this case 
the distribution function is a step function with the "^r-onic fermi energy" 
/^q, Ug^ = gq~^0{iiq — {q — h, s) where = uo^ and = e^. The energy 
is 

3 7 

E = -Nd(2tXsXh + -j^l) (3.1) 

where the order parameters, Xq {q = ^5^), are implicit functions of the 
statistics, a. In Figs. 4 (a) and (6), the energy per lattice site as a function of 
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a for several values of doping rate are shown for Casel and Case2. One can 
easily see that for both cases slave fermion (a = 1) is favored for low doping 
region and slave boson {a = 0) for high doping region. However there exists 
a difference. In Casel the optimum statistics changes continuously from 
slave fermion to slave boson through the intermediate statistics (0 < < 1, 
to be referred to as slave g-on in the following) for 0.351 < 5 < 0.392, as 
seen in Fig. 5(a). On each boundary of two different statistics the transition 
is of second order. On the other hand, in Case2, the optimum statistics 
changes from a — 1 to a = discontinuously at 6 — 0.371, a first order 
transition as seen in Fig. 5(b). This first order transition is accompanied by 
a discontinuous change of the fermi surfaces of holons and spinous. 

The result can be understood as due to bose-condensation of bosons. 
In the low doping region, spinous are abundant and govern the energetics; 
the energy gain is maximum if all spinous are bosons and condense into the 
bottom of the band. Hence, slave fermion (o; = 1) is favored at low doping 
as also has been shown in Ref. 33. On the other hand, in high doping region 
holons are abundant and the maximum energy gain occurs if they become 
bosons and condense and slave boson {a = 0) is favored. 

3.2) Finite Temperature 

Let us first study high temperature behaviours. The onset temperature 
of the w-RVB state is 

T. = ^.„4l + ^l + 2^(§)^], (3.2) 

where aq = 6q{l - 6qgq){l + 6qgq), gq = 1 - gq , {q = s), 6s = ^(1 - 6) 
and Sh ^ S . This expression is derived by assuming that the transition is 
of second order, and does not hold when the transition is of first order (see 
§3.3, and §3.4). 

Above this temperature the effective hopping integral or bandwidth is 
zero and there is no energy gain due to itineracy. Then, the free energy for 
T > Td is determined purely by the entropy, 

S = N{6ln6 + (1 - 6gh)Hl - dgn) - (1 + (5^7.)ln(l + (5^7.) 
+2[5s\n5s + (1 - 6sgs)H'^ - SsQs) - (1 + SsQSM'^ + SsgS)]}. (3.3) 
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Fig. 6(a) and (b) show —S/N = F/NT (free energy per lattice site divided 
by temperature) for different values of doping rate for Casel and Case2, 
respectively. These fix the feature of the phase diagram at T > To; the 
intermediate statistics region is possible for 0.371 < 5 < 0.457 in Casel 
but not possible in Case2. We remark that this is independent of spatial 
dimensionality, since Eq.(3.3) is. As the temperature is lowered through 
To, the transition to li-RVB state is of second order for both cases, except 
for a small region of doping 0.4043 < S < 0.4124 (from A to B of Fig. 7(b)) 
in Case2 where the transition is of first order. This is because the Td's on 
both sides of slave fermion-slave boson interface do not fall on the same 6, 
i.e. the Td from A to B coincides with the interface of slave fermion and 
slave boson [see inset of Fig. 7(b)]. The overall phase diagram for Casel and 
Case2 are shown in Fig. 7(a) and (b), respectively. 

In Casel, the width and the location (in 5-space) of the region where 
the slave g-on is stabilized, is almost independent of T because —S/N, as 
a function of a, has basically the same features as E/N for the same range 
of doping S. 

In general, the absence or presence of intermediate statistics can be 
understood in terms of the total derivative of the free energy with respect 
to a, 

dF{a) 



da 



= T[g'{a)fg,-2g'{l-a)fg^] (3.4) 



where fg = ^gnc,ln(l + Wg ^), {g = gh,9s) and the prime denotes differ- 
entiation with respect to its argument. For Casel, at low doping ^^^^ is 
negative while at high doping it is positive for any value of a. This is due 
to fact fg^ > fg^ {fgy^ > fgj foY low (high) doplug. At the same time it is 

seen that < and > for some range of 5. With this we note 

that there exist a region of 6 where intermediate statistics is possible. 
For Case2, on the other hand, it is straightforward to see that 

has only one zero for < a < 1 and ^^^^ > and ^^^^ < 0, thus, 
intermediate statistics is impossible to occur in this case. 

We also considered another hypothetical assignment of statistics g = 
o? where at — and a = \ reduce properly to boson and fermion, 
respectively. In this assignment, a wide region of intermediate statistics 
was found because < and > . 

The possibility of intermediate statistics may generally be understood 
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by an effective exclusion statistics Qef / defined by pe/ f = 9h-^ 9s- If 9ef / > 

1, it is impossible to realize, whereas, if Qeff < 1, the possibility of its 
occurrence is very high. The case of Qeff = 1, which corresponds to Casel, 
is marginal and elaborate analysis is needed. 

3.3) Effect of Dimensionality 

The 9-on as a generalization of boson and fermion is not limited to 
two dimensional in contrast to anyon whose realization is restricted to d = 

2. In this subsection we extend our calculations to Id and 3d under the 
assumption that the free energy given by Eq.(2.9) is still valid in these 
cases and explores its consequences that may be relevant in understanding 
more on the '^f-on physics' for the t-J model. 

In this subsection only Casel will be discussed where the role of di- 
mensionality is profound in contrast with Case2 which is insensitive to the 
change of dimension. The phase diagrams for Id and 3d are shown in 
Fig. 8 (a) and (b), respectively, where the bose condensation temperature, 
Tb, is also drawn for 3d. The main effect of increasing the dimensionality is 
to reduce the slave g-on region as seen in Figs. 8 (a), 7(a), and 8(b) for Id, 2d, 
and 3d, respectively. This is also seen in Fig. 8(d) where, at T = 0, the op- 
timum statistics as a function of S for each dimension is shown. The energy 
as a function of a for different values of 6 in each dimension further sup- 
ports these understanding as seen in Figs. 9(a) (Id), 4(a) (2(i) and 9{h){3d]. 
Clearly, this suppresion of slave ^f-on-phase as the dimension increases is 
due to increasing tendency of bose condensation. 

In Id, at T = 0, the entire region of 6 favors slave ^-on-phase except 
at the endpoints, 5 = and 6 = 1 where slave fermion and slave boson, 
respectively, are realized. On both ends the transitions are of second order. 
However, the deviation of optimum statistics from a = 1 near 6 = and 
from a = near 6 = 1 are very small as seen in Fig. 8(c). At finite temper- 
ature, in both regions of doping rate the change of statistics are continuous 
also. However, even at very low temperature the region of ^r-on is reduced 
appreciably compared to that at T = as seen in Fig. 8 (a). This is due to 
singularity of the density of state at the band bottom which make the effect 
of entropy important even at very low temperature. 

In Sd, on the other hand slave boson and slave fermion phases are 
stabilized much compared to 2d. This is due to bose condensation. The 
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slave g-on phase is only realized at temperature near and above To where 
the effect of entropy is significant. We note there exists a triple point 
Ptr(^tr — 1.732J, (5tr-=0.373) where the three phases (slave fermion, slave 
^r-on, and slave boson) coexist, at relatively high temperature [Fig. 8(b)]. 
This is because Te scales with the band width and has the same order of 
magnitude as To. Below Ptr, the transition at fixed temperature from slave 
fermion to slave boson is of first order. At Ptr, the slave fermion-slave boson 
interface branches into slave fermion-slave ^f-on and slave ^r-on-slave boson 
interfaces. It turns out that the transition from slave fermion to slave g-on 
phases (line PcrQ) is of second order, while the transition between the slave 
g-on and slave boson is of first order on the line PtrPcr but becomes of the 
second order on Pcrf^- In Fig. 8(c), the a-dependece of the free energy are 
shown at X, Y, Z (inset) where X and Y are points in slave fermion-slave 
boson boundary and Z is inside the slave g-on region which is chosen such 
that the free energies of slave fermion and slave boson are the same. 

Another consequence of Tb<Td is that bose condensation of spinous 
becomes of first order for 6 < 0.0530. This region is indicated by the dotted 
line in Fig. 8(b). Moreover, we found numerically that for very small S (< 
0.0062), the transition into the li-RVB state occurs at higher temperature 
than that given by Eq.(3.2), which has been derived by assuming second- 
order transition. What actually happens is that 1/--RVB is driven by bose 
condensation of spinous, namely. To coincides with Tg, and becomes of first 
order. This is seen from lines 1 and 2 in Fig. 10(a), where the Ginzburg- 
Landau free energy, SF = F{xs)—F{xs = 0), as a function of Xs has two 
minima. Note that the bose condensation is stabilized for Xs larger than 
the value indicated by dot in each line. At slightly larger 5(>0.0062), the 
w-RVB transition becomes of second order but bose condensation remains a 
first order transition (lines 3-7 in Fig. 10(a)), until 6 = 0.0530 beyond which 
Tb becomes of second order. The phase diagram very close to half filling is 
shown in Fig. 10(b). 

Hence, as an overall view the richness of our phase diagrams discloses 
the intricate interplay between bose condensation, band edge singularity of 
the density of states, and entropy at high temperature. In Id, entropy is 
the key factor, except at T = 0. However, in Sd bose condensation is the 
key factor. In 2d, the intermediate statistics are realized on the same range 
of both high and low doping and then we may infer that the tendency to 
bose condense and the effect of entropy have the same weight. 
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3.4) Effect of J 



We next study the effect of J. The global features of our results do not 
depend on the magnitude of J . For example, the phase diagram for J = 
in 2d is shown by the dotted lines in Fig. 7 and do not deviate much from 
the results of finite values of J. The same is true for Id and Sd. However 
in 3d for J = 0, the transitions at Td and Tb are both of second order at 
low doping in contrast to the case of finite J. Moreover, in Casel (Case2) 
the Tb starting from the interface of slave 5f-on(slave fermion) -slave boson 
phases up to ^ = 0.6439 is of first order transition. Qualitatively, the main 
effects of J > are pronounced at low doping region where the number 
of spinous are large. For example at half-filling, it removes the degenaracy 
in the optimum statistics at T = and supports finite To, as well as Tb 
in 3d. Furthermore, it stabilizes the slave fermion-phase and, in general, 
raises T^ especially at low doping region. In Casel, the former results in 
the narrowing of the region of the intermediate-statistics phase. All these 
can be understood as the result of increase in the coherency of spinous, and 
consequently of holons, due to J. 

4. Summary and Discussions 

In the slave-boson theory of the t-J model, the electron operator is de- 
composed into spinous and holons, which are assumed to be fermions and 
bosons, respectively. Within the mean field approximation, the spinous and 
holons are assumed to move independently and the characteristic phase di- 
agram is predicted as given by Fig.l. The most important fluctuations 
around this mean-field solution are the gauge fields, which by themselves 
couple to spinous and holons leading to the effective interaction between 
them. Especially once the bose-condensation of holons is realized, i.e. 
T < Tb, these effects of gauge fields are essential resulting in the bind- 
ing of spinon and holon which implies that the ordinary view based on 
electrons is recovered. In the region where the spin and charge are sep- 
arated, i.e. in the region of anomalous metal (AM) and spin gap (SG) 
in Fig.l, various predictions have been made for spin excitations based on 
the mean-field approximation without the explicit recourse to the fluctu- 
ation effects, which are expected to be relatively minor because of large 
characteristic energy due to Fermi degeneracy of spinous. Above all the 
possibility has been pursued® '^^ of the identiflcation of the spin gap phe- 
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nomenon with the onset of the SG phase out of the anomalous metaUic 
(AM) state, the former and the latter are singlet and uniform RVB states, 
respectively. This identification has been first suggested by Rice^'^^ who 
payed attention to the quantitatively different temperature dependences of 
the NMR shift and rate in the optimal and underdoped YBCO, and this 
fact has been further explored in Ref. 35. This identification has led to 
the prediction of the anomalous frequency shift of the particular in-plane 
oxygen phonon modes^^ , which has been confirmed experimentally^^^ . The 
results of the recent angle-resolved phtoemission spectroscopy seem to be 
also in accordance with the prediction^^\ 

In contrast to these spin excitations, the charge transport properties, 
i.e. the resistivity and the Hall effect coefficient, are very susceptible to the 
fluctuations around the mean field solutions, since the holons are assumed 
to be pure boson. Moreover there are theoretical indications that these 
fluctuations can be very large. 

The present study is intended to overcome this apparent difficulty by 
proposing a new theoretical framework to treat the large gauge fluctuations. 
By attaching the gauge fluxes to spinous and holons, which are assumed to 
have a flxed value but with the fluctuating signs (like Ising model), and then 
treating these particles as free ^f-on's, the statistics has been optimized by 
minimizing the free energy on the plane of doping rate and the temperature 
by use of the distribution functions of gf-on's. (The brief report of the present 
study has been given in Ref. 23.) To best of our knowledge, this is the flrst 
theory to search for statistics variationally and moreover this theory repects 
the time-reversal symmetry. The result of the optimization is that the slave 
fermion, where spinous are bosons, is favored at low doping (near half- 
fllling) while slave boson, where spinous are fermions, is favored at high 
doping, with the intermediate statistics in between. This general feature of 
slave-fermions and slave-bosons are easily understood by the existence of 
the bose condensation of spinous at low doping and holons at high doping. 
This understanding also clearly indicates that the regions of pure slave- 
fermions or slave-bosons are overemphasized in the present approximation 
because the existence of hard cores of spinous and holons is totally ignored. 
Consequently we expect the region of the intermediate statistics will be 
enlarged in a more elaborate treatment. In such a region of intermediate 
statistics the gauge fluctuations are greatly reduced because of the existence 
of the finite mass in the gauge propagators due to the Chern-Simons term^^^ . 
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It is possible that such a region of intermediate statistics (but close to the 
slave-boson) is the one where many interesting and anomalous properties 
have been observed experimentally. 

At the same time we also note the following questions which still re- 
main to be addressed properly. The assignment of statistics, g = g{a), 
we adopted in the present study is derived for anyons in continuum space 
and does not reflect the actual dispersion of our spinous and holons. More- 
over, the assignment was been made by comparing the second virial co- 
efficients which are appropriate only at high temperature. Hence, further 
investigations are needed on the lattice effects as well as the temperature 
dependences. 
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Note added: After submitting this paper, we found that Huang^^^ gave an 
another description of ^'-on in terms of ensemble of bosons whose fraction g 
is randomly transmuted to fermions. Based on his description an alternative 
formulation of our theory is possible; gauge flux absorbed by spinous and 
holons is fluctuating discretely with values, and tt,"^^^ and thus their local 
statistics is fluctuating between slave boson and slave fermion. Note this 
picture applies to any dimensions. This description leads to the assignment 
of statistics deflned as Casel in the text. 

Figure Captions 

Fig. 1. The slave boson mean-field phase diagram. The phases realized are 
antiferromagnetic (AF) state, superconducting (SC) state, "spin gap" 
(SO) phase, anomalous metallic (AM) state, and electron liquid (EL) 
state. 

Fig. 2. Temporal view of spatial variation of gauge fiux (a), its attachment to 
particles (b), exchange statistics (c) and exclusion statistics (d). 

Fig. 3. Distribution function of g-on, n{£), where e, /i and f3 are energy, chem- 
ical potential and inverse temperature, respectively. 
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4. The energy as a function of a for several values of 5 for (a) Casel and 
(b) Case2. In Casel, the inset shows those at the other values of 5. 

5. The optimum values of a at T = as a function of 6 for (a) Casel and 
(b) Case2. 

6. The entropy for T > as a function of a for a choice of several 
values of 5 for (a) Casel and (b) Case2. These are independent of 
temperature, spatial dimensionality and J. 

7. The phase diagrams in d = 2 for (a) Casel and (b) Case2. The inset 
in (b) shows the region where the transition to li-RVB state is of first 
order (from A{C) to B{D) for J = t/4 (J = 0)). The solid and dashed 
lines are for J = t/4 and J = 0, respectively. 

8. The phase diagrams for (a) Id and (b) 3d. In (b), the dotted line 
in Tb and in Tjj represents first order transitions, (c) The enlarged 
version of (b) where slave g-on phase is realized. The inset shows the 
a-dependence of the free energy measured relative to slave fermion or 
slave boson on points X, Y, Z. (d) The doping dependences of the 
optimum statistics are shown for d = 1, 2, 3 at T = 0. 

9. The Q!-dependence of energy for several values of 5 in (a) Id and (b) 
3d at r = 0. 

10. (a) The Ginzburg-Landau free energy as a function of w-RVB order 
parameter, Xs, close to half-filling at temperature chosen such that the 
two minima have the same value. At the dot mark on each line bose 
condensation sets in. Lines 1-2 show why To coincides with Tb and 
becomes of first order transition. The series of lines 1-7 shows how the 
two minima of Xs smoothly merges into single point, i.e. Tb reduces of 
second order as 5 increases, (b) The 5-T phase diagram or small values 
of 5 where the Td coincides with Tb. The transition in Td changes 
from first order to second order t 6 = 0.0062. 
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